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. It has been proposed that propagation of cosmic rays at extreme-energy may 
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be sensitive to Lorentz-violating metric fluctuations ("foam"). We investigate the 
changes in interaction thresholds for cosmic-rays and gamma-rays interacting on 
the CMB and IR backgrounds, for a class of stochastic models of spacetime foam. 
The strength of the foam is characterized by the factor (E/Mp)"", where a is a 
^ I phenomenological suppression parameter. We find that there exists a critical value 

^ ! 

CN ■ of a (dependent on the particular reaction: acrit ~ 3 for cosmic-rays, ~ 1 for gamma- 

. rays), below which the threshold energy can only be lowered, and above which the 
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threshold energy may be raised, but at most by a factor of two. Thus, it does not 



' appear possible in this class of models to extend cosmic-ray spectra significantly 

^ • beyond their classical absorption energies. However, the lower thresholds resulting 
<D '. 

r~| , from foam may have signatures in the cosmic-ray spectrum. In the context of this 



foam model, we find that cosmic-ray energies cannot exceed the fundamental Planck 
scale, and so set a lower bound of 10^ TeV for the scale of gravity. We also find 
that suppression of p ^ pn^ and 7 — > e~e+ "decays" favors values a ^ Ocrit- Finally, 
we comment on the apparent non-conservation of particle energy-momentum, and 
speculate on its re-emergence as dark energy in the foamy vacuum. 
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I. INTRODUCTION: METRIC FLUCTUATIONS AS SPACE-TIME FOAM 

The space-time metric tensor g^i, becomes a dynamical variable when gravity is quan- 
tized, and space-time foam is the quantum mechanical uncertainty 5gfj,i, in this variable. We 
investigate the possibility that the foam has phenomenological consequences. 

We work within the framework of foam models having a parameterization given by 



where Zp = y ^ is the Planck length and tp = ^ is the Planck time in a four-dimensional 
theory. In theories with more than four dimensions, these scales could be larger, in fact, 
much larger, than the Planck scales. The parameter "a" depends on the foam model [1, 2]. 
The covariance implicit in the fluctuation variable Sg^i, puts space and time uncertainties on 
an equal basis, in contrast to the situation in non-relativistic quantum mechanics where the 
length-momentum relation arises from operator commutators, and the energy-time relation 
from a less-compelling argument using the Fourier decomposition. The uncertainty relations 
given in Eq. (1), and those that follow below, are sometimes called "Generahzed Uncertainty 
Principles" (GUP). We will use this name. 

From (1) and the fact that Sf = l^Sg, one obtains the uncertainties for length and time: 



The quantum mechanical relations between length and momentum, and time and energy, 
then lead to the equivalent expressions for the uncertainty: 



where Mp ~ lO^^eV denotes Planck mass. From here on we drop explicit mention of powers 
of c. These uncertain lengths, energies, etc., make the lengths of four- vectors uncertain, 
and so break Special Relativity. Accordingly, we must single out a frame in which these 
uncertainties are defined. It is common to assume that the special frame is the cosmic 
rest frame, in which the cosmic microwave background (GMB) is isotropic. We adopt this 
assumption. 




(1) 







(3) 
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The uncertainty relations defined above can effect violations of Lorentz invariance (LIV) 
even in the weak-field, flat-space hmit. Two of the most studied cases are (i) modification of 
the energy-momentum conservation equations, AP'^ = 0, or (ii) modification of the energy- 
momentum dispersion relation, = irP'. Of course, LIV may also appear in both (i) and 
(ii) simultaneously. In this work, we study LIV of the energy-momentum conservation law 
at the scales defining the GUP, but maintain the usual dispersion relation — pj + rnf. 

II. MODIFIED THRESHOLD ENERGY FOR UHECR ABSORPTION 

Following many others [1]-[10], we investigate the role that LIV-kinematics may play on 
ultra-high energy cosmic rays (UHECRs). The UHECRs are gamma-rays at ~ tens of TeV, 
and nucleons at E 10^° eV. Standard particle and astrophysical arguments lead one to 
expect TeV gamma-rays and 10^° eV nucleons to be just above their respective thresholds for 
absorption on cosmic radiation background fields. The annihilation reaction for TeV gamma- 
rays is 7 + 77i?B— >e"'"-|-e~, where ^ymB denotes a cosmic infrared background photon; the cms 
energy threshold is ^/s — 2me. The energy-loss reaction for 10^° eV nucleons is dominated 
by + 7cMB— N' + 7: near threshold, where 7cms denotes a photon in the CMB; the 
cms threshold energy is ^/s = m^. 

The generalized uncertainties can in principle raise or lower the energy thresholds of 
these reactions. Raising an annihilation or energy-loss threshold presents the possibility of 
extending the CR spectrum beyond expected cutoff energies. There have been suggestions 
that such extended spectra do exist for both gamma-rays [8] and for nucleons [9] . Predictably, 
there have also been suggestions that LIV is the origin of the anomalous spectra [2] . 

Consider the general 2 —> 2 scattering reaction a + 6— >c + d. Let ma, mi,, mc and m^ label 
the particle masses, and Ea, Ef,, Ec and E^ the particle energies. The unmodified dispersion 
relation reads 



However, the energy and momentum conservation laws including GUP fluctuations become 





Ea + SEa + Eb = E, + 6E, + Ed + 6Ed , 
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Pa + ^Pa +Pb ^ Pc + 5Pc + Pd + ■ (5) 

According to Eq. (3), the energy and momentum fluctuations of the cosmic background 
photons, 5Ei, and 5pij, are very smaU since and are so smaU (~ meV for the CMB and 
a few tens of meV for the far-IRB). We have set them to zero. The magnitudes of the other 
uncertainties can be significant. 

Subtracting the second of Eqs. (5) from the first, inserting (4) for the differences E — p, 
and reafizing ^ that Ph = —E^, one gets for the energy of the background photon, 

= El + ^Foam , (6) 

where 

^\Pc Pd Pa] 

is the classical value, and 

(^Foam = 4 + (^d - (^a , with 8j = ^ {5Ej - 5pj) , (8) 

is the contribution from fluctuating foam. To flrst order in m|, the pj's in Eq. (7) can be 
replaced by Ej^s. We remark at this point that (^poam is sourced by SE and/or Sp. Therefore, 
the origin of ^poam is open to a broad interpretation. We discuss this a bit more in section 
(VI). Also, since 6p is a component of a three- vector, whereas 6E is not, the relative sign in 
the combination 5E — Sp appearing in the deflnition of ^Poam is not meaningful. 

Next, consider the reaction at threshold, ^/s = rric + rrid. In terms of the boost factor 
7 between the center of mass frame and the lab frame, one has E^^^ — j{mc + m^) = 
Ea + 0{Eb), Ec = -frric, and E^ = -fma, i.e., 7 = |^ = |^ = T^^mi' ^^^^^ equalities allow 
the elimination of Ec and E^ in Eq. (7) in terms of Ea, which we write as E^^ to remind 
ourselves that the kinematics are being calculated at threshold energies. The result is that 
Eq. (6) becomes 

4 Eb -Eth = {rric + m^)^ -^1 + E^y, ^poam • (9) 

^ we choose to work in a frame where CR momentum is pcR = Pcr{Q, 0, z), and momentum of background 
photon is pb = Pb{0, 0, —z) 
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Solving this equation for E^^ then gives the modified threshold energy for the reaction. Of 
course, some model for the fiuctuations must be introduced. 

References [1-3] argued for the following form of the fluctuations: 

5,^l-iSE,-5p,)^-'-p, (^Y^~E, (^y. (10) 



Different choices for a and e parametrize different space-time foam models. With dimen- 
sionful mass-energy factors explicitly shown, e is expected to be a number roughly of order 
one. The exponent a is assumed to be positive such that fluctuations are suppressed below 
the Planck scale. Possibilities other than the particular form for the fluctuations given in 
Eq. (3) are certainly possible. For example, one may choose to parameterize the fluctuation 
(^Foam in Eq. (9) directly, with a parameterization of one's choosing. However, we will stay 
with the form given above. 

As we have seen, the individual Ej^s are linearly related to each other by mass ratios. 
Along with Eq. (10), this means that the S/s are also (nonlinearly) related to each other by 
mass ratios. This makes the result of inserting Eqs. (10) for each j = a, ,b ,c into Eq. (9) 
fairly simple. The result is a general and manageable equation for the reaction threshold 
energy, incorporating the correction from space-time foam: 

n2 f Eth 



4 Eb -Eth = 4 £'b -Eciass + e -E", 



th 
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with -Eciass = (™'^+^^|^ being the energy of the threshold when special relativity is not 
violated. When e = 0, the classical threshold £^th — -Eciass of course obtains. When a — 0, 
the classical threshold E'th = -Eciass also obtains, for any value of e. 

Before examining speciflc models for the meaning of a and e, we can extract from Eq. (11) 
the number of real positive roots of E^^. These roots are the candidate solutions for the 
modified threshold. Imagine plotting the LHS and RHS of the equation versus -Eth- The 
LHS of (11) rises linearly in E^^ from zero, with a slope of 4Eij. The RHS rises (falls) for 
positive (negative) e at a higher power of E^th, from a positive intercept of {rric + rudf' — m^. 
For negative e, the two curves will always cross once and only once, i.e., there is always a 
single positive root. For positive e, the two curves may never cross, "kiss" once, or cross 
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FIG. 1: Schematic illustration of the classical and modified threshold energies. 

twice, giving none, one and two positive roots, respectively. There is a critical value of 
a = Ocrit, dependent on the particle masses and the positive value of e, at which there is a 
single positive root, above which there are two, and below which there are none (since a is 
the exponent of a ratio less than one). These results are illustrated in Fig. (1). 

When e < 0, the single positive solution for Eth is lower than the classical value -Bdass- 
This leads to cutoffs in the CR spectra at energies lower than those predicted from classical 
physics. When e > and a < Ocrit, there is no physical solution for i?th, and so the absorption 
reaction does not happen at any energy. We return to this case briefly in §V. When e > 
and a > Ocrit, then the solutions for £Jth are always larger than £^ciass- This leads to CR 
cutoff energies higher than those predicted from classical physics. However, as a increases, 
the influence of the foam term decreases, and for a > Ccrit, the lower of the two solutions, 
call it E^^, approaches the classical value £^ciass- The higher of the two solutions, call it E^, 
goes to Mp as a increases. 



III. FOAM DYNAMICS 



Not surprisingly, there exist several models for the foam dynamics of a and e. We discuss 
some of these in this section. 
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A. 



Foam Model with Fixed Fluctuation Parameter e 



Since one interest of the particle-astrophysics community is to explain possibly extended 
CR spectra, we first discuss the higher cutoffs provided by the Eth > -Bciass case, i.e., the 
e > 0, a > If there were a reason in Nature to reject the lower E^^ solution (we 

know of none), then the arbitrarily-large value of Eth = E^ would allow CR cutoffs to be 
arbitrarily extended. However, with both solutions operative, the reaction will occur when 
E rises to exceed either solution, i.e. to exceed mm{E^, E^^^}, which is just the lower-energy 
solution E^; the higher E:^ solution seems irrelevant. So, how large can E^ be? The answer 
is that E^ is maximized at the single solution value occurring when a — Ocrit- Call this value 
^max These solutions and their labeling are shown in Fig. (1). 

The critical a and the critical Eth can be found in principle by simultaneously solving 
two equations.^ The first is just (11), and the second is obtained from (11) by equating the 
first derivatives of the LHS and RHS with respect to Eth- However, manipulation of these 
two equations does not lead to a useful analj^ic separation of Ocrit and E^^. We content 
ourselves to use numerical techniques in the main text to determine Ocrit and E!^^, but 
present some accurate analytical approximations in the Appendix. There is however, one 
simple analytic relation that results from manipulations of these two equations. It is 



This result shows that E^^ depends on e (assumed positive here) only implicitly through 
Ocrit, and that E^^ lies in the interval [E'dass, 2 E'dass] regardless of the value of Ocrit- 
approaches 2 E'dass as a^rit — ^ 0"*", and approaches ^ciass for acrit 3> 2. There are E^ solutions 
exceeding 2 Edass, but these are inevitably accompanied by a second solution, E~^^, lying 
below 2 Eth- This is our first new result. We repeat it: For positive e, the reaction threshold 
energy can be raised, but at most by a factor of 2. 

Let us comment on the foam- inspired extended CR spectra obtained in [1-3]. In this work 




(12) 



^ In principle, Eq. (11) is sufficiently nonlinear that there could be more than one solution for the critical set 
{Eth, a}- However, more than a single solution is not apparent in our numerical work. In our Appendix 
we explain why the single solution for acrit is effectively unique. 
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a pre-desired value of E'th is input into Eq. (11), and from this the value of fixed, positive 
e is extracted. This approach suffers from (at least) two drawbacks. The first is that it is 
oblivious to the existence of the second, lower-energy solution E^. The second drawback is 

that it is fine-tuned in the value of e. The first drawback is much more serious, for we have 
just shown that the second solution raises the threshold energy by at most a factor of 2. The 
model of [1-3] appears to fail to raise threshold energies. 

B. Foam Model with Gaussian Fluctuations 

A healthier approach to foam dynamics is described in [4]. The parameter e is treated 
as a stochastic variable, subject to some specified probability distribution.^ The stochastic 
assumption seems reasonable, in that space-time fiuctuations at one site would not depend 
on the fluctuations elsewhere. We follow this approach here, generalizing the results of [4] . 

When an interaction occurs, the kinematics are determined by the fixed parameter a and 
a single random value of e. Since experiments sum over many events, the total data sample is 
best described by the most probable value of the threshold energy. This is determined by the 
Gaussian-distributed e. In principle, some random occurrences will reduce the threshold for 
particular events even below the mean threshold. However, our numerical work reveals that 
the width in i?th, resulting from the distribution in e, is small. This is evident in Figs. (4)-(7). 

In principle, the value of the fiuctuation Si of each particle can be treated as independent 
stochastic variables, the effect of this on Eq. (11) is to replace the overall e with an inde- 
pendent Ea, ec, ed for the the respective three terms in the bracket. In the end, an overall 
Gaussian is probably a good approximation even for this more complicated case, given the 
generality of the central limit theorem. One might expect quantitative, but not qualita- 
tive, differences [4]. We will treat the overall e as the single stochastic variable, and define 
the probability distribution of the single e as p(e). We will follow [4] and assume Gaussian 

^ The model of Ng ct al. can be thought of as the special case where the stochastic distribution of e's 
is sharply peaked about a pre-determined positive value and the lower of the two energy-thresholds is 
ignored. 
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statistics. Then 

p(6) = ^e-^. (13) 
(Tv27r 

The stochastic e's are then generated numerically via 

€ = aV2erf-\r) + e , (14) 

where r is a random number in the interval [—1, 1], a is the variance of a distribution, and e 
is the average value of e. We set cr = a to avoid introducing a new parameter. Alternatively, 
one could for example choose a constant variance, say cr = 1, or any other value. This does 
not change the general behavior of our results, as we show in the next section. We choose 
e = based on a preference for symmetry, and to maintain the smallness of fluctuations. The 
same choice was made in [4]. A model with nonzero e was proposed in [5]. Obviously, this 
expresses a preference for negative e (lowered threshold) over positive e (raised threshold) or 
vice versa. While this asymmetrical choice may turn out to merit Nature's attention, it has 
not yet attracted our attention. 

With the symmetrical choice for e, half of the fluctuations present negative e, and half 
present positive. For the negative half, each e generates one solution for Eth, with E'th < 
-E'ciass- For the positive half, each e generates no solution when a < Ocrit, and two solutions 
above fi'ciass when a > Ocrit- The lower of these two solutions, E^, is relevant, while the 
higher solution is probably not. 



IV. MODIFIED THRESHOLDS IN DETAIL 



For gamma-rays incident on the IRB, Eq. (11) becomes 



(15) 



where rrie is the electron mass, and for definiteness we take Eirb — 0.025 eV. For CR 
nucleons interacting on the CMB, Eq. (11) becomes 



AEcMBEti, = {rrip + m^f -ml + e 



(mp + m^)i+" 



(16) 
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FIG. 2: Most probable values of threshold energies (in eV) vs. foam models for 7-rays 

with nip and the nucleon and pion masses. Here, for definiteness we take Ecmb — 
7.2 X 10~^ eV, near the mean energy of the spectrum. Solving these equations numerically, we 
map the random Gaussian distribution of e described in Eq. (14) onto a random fluctuation 
spectrum for i?th- 

Different choices of a characterize different foam models. The two choices a = | and 1 
are motivated by interesting plausibility arguments [2, 3]. Integral values of a are motivated 
by loop quantum gravity [6], and also by arguments for unbroken rotational invariance [7]. 
We take the agnostic approach and treat a as a continuous parameter to be explored from 
zero upward. In Figs. (2) and (3) we show the evolution of solutions with a. Each "solution" 
£^th is really the most probable value of Et\i picked from a distribution.^ 

The values of a^^t are evident in the critical points of these figures. Numerically, they are 
Ocrit — 0.964 for the gamma reaction, and acrit — 2.87 for the nucleon reaction. Numerical 
values of Ocrit in general depend on the variance a in the p{e) distribution, here set for 
simplicity to a. When cr is taken as a free parameter, we find that the following limiting 
values for Ccrit would result: Ccrit — when cr — > 00, and Ccrit — 00 when o" — > 0. 

Below Ocrit, the single curve reveals the single solution accompanying negative e fiuctua- 

* Numerical work reveals that for e < (e > 0), the most-probable value for iJth is given by the solution to 
Eqs. (15), (16) when e is set roughly equal to —a for e < 0, and set to +(j for e > 0. This is sensible, as 
restricting the two sided distribution p{e) to just one side moves mean e from zero to iv^^TTcr. 
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FIG. 3: Most probable values of threshold energies (in eV) vs. foam models for UHECRs. 
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FIG. 4: Distributions of threshold solutions for gamma-rays 

tions. Above Ocrit, and accepting both positive and negative e fluctuations, three solutions 
are actually present. The two relevant solutions, Et^ from negative e fluctuations, and 
from positive e fluctuations, are nearly identical in value, clustered just below and just above, 
respectively, the classical solution £'ciass- In the figure, these two near-classical solutions con- 
stitute the unresolvable horizontal branch of the curve to the right of the critical point, while 
the irrelevant solution E^ is the curve that rises to the right acrit- 

In Figs. (4) and (5), the individual solution "packets" arc shown, for different values of 
the parameter a. The area of each packet reflects how often a fluctuation in e produces 
a physical solution. Below Ccrit, the total area in the packet is 50%, since only half of the 
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FIG. 5: Distributions of threshold solutions for UHECRs 



fluctuations, the e < ones, produce a physical solution. At acrit, the sharp packet has total 
area of 100%, reflecting one solution for each e < and one for each e > 0. Above Ocrit, 100% 
of the area remains in the sharp peak labeled a > Ccrit, coming from the e < solution and 
the solution, plus another 50% area exists in the higher-energy solution. Raising a 

reduces the magnitude of the space-time fluctuations, and so pushes the packets comprised 
of the e < and E^ solutions closer to £^ciass- Raising a also pushes the E^ solution ever 
higher, toward Mp. 

We note that the width of the e < packet decreases as a — > a~;^. Above Ccrit, the 
E^ width and e < solution width remain narrow, the former related to the bounding of 
E^ between £^ciass and 2£^ciass- The E^ width above Ocrit gets narrower as a increases, or 
equivalently, as E^^ Mp, due to the suppression of the fluctuation. 

In Figs. (6) and (7) we show the Eth packets for gamma-ray reactions and nucleon reac- 
tions, for a = I and a — 1. These two values are among the most popular in the literature. 



Note that a = | is below Ccrit for both nucleon and gamma-ray reactions. Accordingly, 
there is a single solution for Et^ from e < and none from e > 0. The e < solution lowers 

-Eth below the classical threshold. On the other hand, the a = 1 value is again below Ocrit for 
the nucleon reaction, but above acrit for the gamma-ray reaction. Thus, the a = 1 gamma- 
ray distribution shows the nearly-classical threshold, as well as the irrelevant higher-energy 
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FIG. 6: Threshold Distributions for 7 rays (energy in eV) 
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FIG. 7: Threshold Distributions for UHECRs (energy in eV) 
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V. AN ASYMMETRIC FOAM MODEL THAT DOES EXTEND THE SPECTRA 

We have shown that one can only lower the reaction thresholds when fluctuations to 
negative e are allowed. We have also shown that positive e fluctuations raise the threshold 
energy by at most a factor of two when a > Ocrit, but yield no physical solution for E'th when 
a < Ocrit- Thus, there is but a single way to raise the reaction thresholds: the stochastic 
variable e must be restricted to positive values, and a must be below Ocrit; the absence of a 
physical threshold in this case means that the absorption reaction cannot happen at all. 
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To postulate that a < Ccrit seems acceptable. After all, the values for a suggested by 
modelers are small, as we have discussed. The concomitant requirement that stochastic 
e > 0, i.e. p{e) = for all e < 0, seems harder to justify. A Gaussian distribution cannot 
deliver this, but any distribution of finite extent, e.g., a simple top-hat distribution or S- 
function distribution, can. Perhaps Nature is this asymmetric, but single-sided distributions 
for fluctuations seem to us to be contrived, and we honorably choose not to pursue them 
here. 

VI. DISCUSSION, WITH SPECULATIONS 

In this section we present some additional interesting issues of the foam model. 

A. TeV-scale Gravity? 

The numerical results we have presented depend on the assumption that the Planck mass 
has its usual value 1.2 x 10^^ eV. Modern thinking, inspired by the extra dimensions available 
in string theory, admits the possibility that the fundamental mass of gravity may be much less 
than this value. One may ask what changes in our analysis with this (much) smaller mass- 
scale. Qualitatively, nothing changes. The construction of solutions for positive or negative 
e values, illustrated in Fig. 1 is unchanged. Also, the relation between El^^^ and -Eciass given 
in Eq. (12) is independent of Mp, and so is unchanged. So it is just the quantitative values of 
the new threshold energies that are changed. But the quantitative changes can be dramatic. 

In our numerical work we have investigated the effect of lowering Mp. We find that, except 
for very small values of a, E'th cannot exceed Mp. This means that the Planck mass Mp 
provides an energy limit for all cosmic rays. Turning this remark around, the observation of 
cosmic rays above 10^" eV implies, in the context of this LIV model, that Mp > 10^° eV. The 
bound is eight orders of magnitude below the usual 4-dimensional Planck mass. However, 
the bound is also incompatible with the popular TeV-scale gravity models by about eight 
orders of magnitude. 
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B. Proton and Photon "decays" 

The general form of Eq. (11) enables us to discuss other processes forbidden by exact 
Lorentz invariance. Above some threshold energies, certain particle "decays" may become 
kinematically allowed [4]. In [4, 5], the authors discuss the two specific decays p pn^ 
and 7 — > e~e"'" for the fixed value a = 1. In this section we will follow [4] and assume that 
energy-momentum fluctuations of the initial state particle characterize the Lorentz violation. 
However, we generalize their results by treating a as a continuous free parameter. Following 
a treatment similar to that which led to Eq. (11), we here find the energy-threshold for the 
process a — > c -|- d to be 

Etn = { [(me + maf - m^] l^j'^^"^'^ . (17) 



One may invert this equation to isolate the a-parameter: 

(18) 



a — 



log r^[{mc + mdf -m^^ 

th 



log (i^) 



It can be easily seen from the form of (17) that only for e < is the threshold energy E^^ 
positive and real. For e > there is no physical solution. Furthermore, for any negative e 
and any a, there is a unique Etu- We present our results in Fig. (8) for —e = 1. As with the 
scattering processes presented earlier, here too predicted thresholds are bounded from below 
by some simple function of particle masses (given in Eq. (17)), and from above by Mp. For 
the a = 1 case of [4], one obtains Efh ~ l^^^eV for "pion-bremsstrahlung" by a free proton, 
and -Eth ~ for photon decay to e^e~ . To hide the particle decays, one may move 

the threshold energies beyond the highest observed cosmic ray energies (-Eciass) by raising 
the value of a. Doing this, one gets from Eq. (18) a lower limit for allowed values of a. To 
be specific, inputing 10^° eV for cosmic rays, and 10^^ eV for gamma-rays, one obtains for 
€— —1 the lower bounds of 2.82 and 0.93 for the respective a's.^ 



^ We note that the lower limits for a's from particle decays are very similar to the lower limits for a's from 
the related scattering processes, given by the Ocrit values obtained in §IV. 
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FIG. 8: Threshold energies for p pir^ (higher curve) and 7 — > e e"*" (lower curve) decays as a 
function of parameter a. 

C. Cosmic- Ray "Knees" 

We note that although our results suggest failure for the program which attempts to raise 
cosmic-ray absorption thresholds, the lowered threshold may nonetheless have application. 
It is conceivable that structure in the cosmic-ray spectrum, e.g., the first and second "knees", 
and the "ankle" , are results of lowered absorption thresholds. In particular, if the parameter 
a itself were to trace the matter distribution, along the lines mentioned above for inhomo- 
geneous dark energy, then Et^ for galactic cosmic-rays could be much lower than for the 
extra-galactic component. The net result might be spectral breaks at these lower i^th values 
for cosmic-rays contained in our Galaxy, i.e., "knees". The argument disfavoring this remark 
is that Galactic cosmic rays are probably not contained long enough to undergo absorption 
(assuming, as we have here, that the absorption cross-section does not depend on e or a). 

D. A Random Walk Through Foam? 



There appear to us to be conceptual problems associated with the kind of model studied 
here. For example, it is an assumption in the model that energy-momentum is transferred 
between the particle and the foam at the point of particle interaction. However, energy- 
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momentum transfer between a free particle and the foam must be disallowed, for this would 
present a random wall through Planck-sized domains, leading to, e.g., unobserved angular 
deflections of light. Resolution of these kinds of ad hoc rules, e.g., smoothing of the foam, 

must await a true theory of quantum gravity. In advance of a theory of quantum gravity, 
speculations are allowable, and in the next subsection we provide a few. 



In foam models, energy and momentum are not generally conserved in particle inter- 
actions. Non-conservation is suppressed for interactions at terrestrial accelerator energies 
and below, but possibly becomes noticeable for interactions of extreme-energy cosmic-rays 
with cosmic background fields, or with atmospheric nuclei. Obvious questions are how much 
energy-momentum is missing, and where does it go? To the best of our knowledge, these 
questions are not addressed in the literature. 

In the context of the model analyzed in this work, it is simple to estimate the energy- 
momentum loss of the interacting quanta. As shown in Eq. (6), the energy- momentum added 
to the quanta is Spoam, which, obtained from Eq. (9), is 



This result makes it clear that a lowered (raised) threshold, ^Poam < ((^Poam > 0), and some 
particle energy-momentum is lost (gained). 

Since the numerical value of E'th depends on a, Spoam too depends on a. We may see this 
directly by comparing Eqs. (9) and (11) and to deduce that 



This equation also makes it explicit that the sign of Spoam is the same as the sign of e. Thus, 
particle energy-momentum is lost (gained) when e < (e > 0). 

We have seen that when a < Ocrit, there is no solution for E'th for positive e, and a single 
solution, typically well below £^ciass, for negative e. Thus, stochastic fluctuations will on 
average drain energy-momentum from the interacting particles, when a < Ocrit- For a > Ocrit, 



E. Energy-Momentum Non-conservation 




(19) 




(20) 
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Foam model a=^ 
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^Foam [eV] 

FIG. 9: Energy-momentum transferred to vacuum foam, as a function of negative e (the case where 
£?th is lowered), for a = 1. 

E'th is very nearly £'ciass; and there is negligible energy- momentum transfer. An example of 
significant energy- momentum loss is displayed in Fig. (9), where a distribution of ^poam versus 
negative e is shown for cosmic-rays, with the sub-critical value a = 1. The typical missing 
energy-momentum per interaction is ~ 120 eV, consistent with Eq. (19). As a fraction of 
the interaction energy-momentum, this loss is tiny. Prom Eq. (19), when E'th <^ -Eciass one 
has for the fractional loss, ^PoamZ-Eciass = -4.Eb/Eti, ~ 10-^^(Eb/EcMB)(PeV/Eti,), far too 
small to ever be detected by direct cosmic-ray measurements. 

Where the missing energy-momentum goes is a difficult question. Without a theory of 
quantum gravity, one is reduced to speculation. Perhaps the vacuum of the Universe itself 
recoils in the interaction (similar to Bragg diffraction, or the Mossbauer effect).^ In an 
isotropic Universe, the net three-momentum transferred will be zero, after averaging over 
many interactions. So our Universe is not running away from us. However, the net energy 
in the vacuum will grow with each interaction. Perhaps an integration over the interaction 
history provides the observed "dark energy" of the Universe, with equation of state ^ = — 1. 

Or perhaps the missing energy-momentum is not homogeneously distributed. For exam- 



(D 




^ In the related model of [11], energy- momentum is transferred to the extra-dimensional bulk of string 
theory, via "D-particle" dynamics. 
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pie, cosmic ray interactions with galactic matter may deposit missing energy-momentum into 
the vacuum of galaxies. Then one might expect dark energy to exhibit some "clustering" 
around large-scale matter distributions. Such a dark energy over-density might appear to 
be " dark matter" . 

An alternative point of view is that the missing particle energy is simply not missing. 
For example, since ^poam = ^{^E — Sp), perhaps SE — and only momentum is lost in the 
interaction. The formalism in this paper goes through unchanged in such a case. As for 
the momentum, even though momentum is lost from the particles in each interaction, the 
isotropy of the Universe guarantees that after averaging over momenta directions from many 
interactions, no net momentum is gained or lost by the vacuum. Thus, there may be no net 
transfer of energy-momentum to the vacuum. 

Particle physics in the expanding Universe may provide a useful focus for thought. In the 
expanding Universe, the momentum in a co-moving collection of particles is not conserved 
(it red-shifts), even though the local Einstein equations conserve the energy- momentum 
of interacting particles. This global non-conservation of energy is related to the lack of 
time-translation symmetry in the global, expanding Universe. Perhaps in the foam models, 
translation symmetry is sufficiently broken, or the interaction of the particles with the met- 
ric foam is sufficiently nonlocal, to make the apparent energy-momentum non-conservation 
palatable. Even the concept of local versus global becomes confused when the particle in- 
teraction involves the foam. If the foam fluctuation is included in the "local" environment 
of the interaction, then energy-momentum can be said to be locally conserved. However, if 
the foam is counted as part of the global vacuum, then energy-momentum is not conserved 
locally, and may or may not be conserved globally. 

F. Relation to Modified Dispersion Equation Approach 

Finally, we remark on the relation between modified interaction kinematics, as presented 
in this work, and related work on modified particle dispersion-relations. It turns out that 
the same modified threshold equation is obtained in either approach, with one significant 
difference: the relative sign of the stochastic variable e is opposite in the two approaches 
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[1, 2]. If one includes both modified kinematics and modified dispersion relations, then the 
two foamy effects cancel each other/ 

VII. CONCLUSIONS 

We have analyzed a family of Lorentz-violating, space-time foam models. Fluctuation 
amplitudes are assumed to be stochastic, Gaussian-distributed, and suppressed by the Planck 
mass. These models can be tested, in principle, by searching for anomalous propagation of 
high-energy cosmic-rays and gamma-rays. We have derived an equation for a new (modified) 
threshold energy, in the general case of two particle scattering. As relevant examples, we 
examined the modified energy-thresholds for the reactions N + jcmb^^ N' + tt and 
7 + liRB^G^ + affecting propagation of extreme-energy cosmic rays and TeV gamma- 
rays, respectively. Our threshold solutions can be parametrized by a suppression parameter 
a. For a given particle reaction, there exists a critical value Ccrit of this parameter, beyond 
which foam does not alter the standard predictions. Furthermore, for a < Ccrit, foam does 
alter the reaction significantly, but always lowering, never raising, interaction thresholds. 
We found that Ocrit ~ 3 for the A'" -|- '^cmb^^ N' + n reaction, and Ocrit ~ 1 for the 
7 + liRB^^^ + e~ reaction. 

Previously, Aloisio ct al. [4] investigated one model (corresponding to our a = 1 case) 
for its effect on the UHECR absorption reaction on the CMB. They found that the reaction 
threshold is lowered to a most probable value of £^th — 2.5 x 10^^ eV. Since the a — 1 
used in [4] is below Ccrit, their lowered threshold is consistent with our more general results 
which extend their model to arbitrary values of a. Even for a > Ocrit, we found that the 
threshold can be raised only by a factor of at most two (according to Eq.(12)), and even 
then only if fluctuations are asymmetrically distributed about zero. Thus, it appears that 
spacetime foam models of the kind we assessed are unable to extend the spectra of UHECRs 
or gamma-rays beyond their classical absorption thresholds on the CMB and IR background. 

Of course, this assumes the same form for the fluctuations e in both the kinematics and the dispersion 
relation. 
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respectively. 

Tests of foam models by means other than anomalous cosmic-ray propagation abound 
in the literature. These include, but are not limited to, proton and even photon decay, 
neutron stabilization, anomalous UHECR shower development, .... Although a discussion 
of all these possible signatures is beyond the focus of this work, we did investigate the 
energy-thresholds above which kinematically disallowed 1^2 processes become allowed. 
The thresholds for these processes increase monotonically with the a-parameter. For the 
processes N + ^cmb^^ N' + n and 7 -|- T/rb— >e+ -|- e~, the energy-thresholds can be 
pushed beyond the observed end-points of cosmic-ray and gamma-ray spectra with a values 
(remarkably) similar to respective Ocnt values. 

We have succumbed to temptation and presented a few speculations. One is that the 
transfer of energy between interacting particles and foam may contribute to dark energy, or 
to apparent inhomogeneous dark matter. Another is that lowered thresholds in cosmic ray 
absorption may effect the "knees" (increase in the spectral slope) observed in the Galactic 
cosmic-ray spectrum. 

Finally, we noted that the absorption thresholds for cosmic-rays and gamma-rays, in the 
context of the model, cannot exceed the "Planck mass". Thus, the observation of cosmic rays 
with energies at 10^° eV necessitate a fundamental Planck mass in excess of 10® TeV. This 
model is, then, incompatible with TeV-scale gravity models, by eight orders of magnitude. 

Acknowledgments: This work was supported in part by DOE grant DE-FG05-85ER40226. 

APPENDIX: APPROXIMATE ANALYTICAL EXPRESSIONS 

In this appendix we derive approximate analytical expressions for various quantities ob- 
tained in the paper by numerical methods. Only the lowest order approximations will be 

given; it is straightforward to calculate higher order corrections. 
We start by rewriting Eq. (11) in dimensionless form: 

X = 1 + e6c-"-V(a)x"+^ (21) 
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where 

6^^, c^^, (22) 

/(a) = 1 - 7-^^^—. (23) 

For the photon (7) reaction 7 + ^irb^g~^ + e~ and nucleon (n) reaction + ^cmb 
A^' + TT, the parameters are the following: 

fti 1.22 X 10^9, ;y 1.17 X 10l^ = 1, (24) 
6„ 4.24 X 10^°, c„ 1.25 x 10^ ^„ « 0.149. (25) 

We consider three different cases with fixed e: (1) positive e and critical a; (2) positive e 
and non-critical a; (3) negative e. Then we generalize to a distribution in e. Because of the 
plethora of subscripts and superscripts required in this appendix, we find it convenient to 
replace the subscript "crit" with the compact subscript "*" . 



1. e > 0, a = ttcrit = a* 

To determine the critical value of a, we need to solve Eq. (21) and its derivative equation 
simultaneously. These two equations are: 

x, = l + e6c-"*-V(a,X*+l (26) 
1 = e6c-"*-V(a,)(a, + 2)xt'^\ (27) 

Eliminating from Eqs. (26) and (27), we find 



ebc-'^-'fia.) = )7 : (29) 



(g. + 

(a. + 2)'^*+2 



Eq. (29) cannot be solved exactly for a*. It is, however, very easy to find simple and accurate 
approximate solutions. To find these, note first that as a* increases from zero to 00, the right- 
hand side of Eq. (29) monotonically decreases from | to zero as (a*e)~^, and the left-hand 
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In (1 + z) In z 



side first increases from zero to its maximum and then exponentially decreases back to zero. 
Thus, in general, there can be none, one, or two solutions to Eq. (29). For large b and c (as 
in both 7 and n cases) and moderate e, there are two solutions: a*o ^ 1 and a*i ~ 1. 

To find the smaller root a*o, we expand terms in Eq. (29) to the lowest order in a* and 
obtain 

(30) 

Substitution of numerical values gives ea*07 ~ 3.45 x 10~^^ and ea*on ~ 1.90x10"^^. Although 
these small values justify the approximation, they are too small to be physically interesting 
for moderate values of e. Thus we learn that for the relevant situation of large b and c, there 
is effectively a unique value for Ucrit- 

To find the physically interesting Ocrit = we note that for moderate a* and large c, 
the term c~"*~^ is the fastest changing part of Eq. (29). Evaluating both the right-hand 
side of Eq. (29) and the function /(a*) to be of order unity, we arrive at the approximation 
a*i fti a^^i + a^^^ , where 

+ in?' (^1) 



(«2' + i)' 



Substituting numerical values, we find 



4°! fa 0.93 + 0.028 In e, (33) 



2.78 + 0.054 In e. (34) 



Figs. 10 and 11 compare the lowest and the next order approximations for function 
of e, with the corresponding numerical solutions for the 7 and n reactions, respectively. 

2. e > 0, a > acrit = 0*1 

As a function of x, Eq. (21) with a > a*i and positive e has two roots, xi and X2- The 
term ebc'"'^ f {a) in Eq. (21) is of order unity for a = a*i and rapidly decreases for a > a*i. 
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FIG. 10: Critical value clS cl function of e for the reaction 7+77iiB— ^-e"*" + e . The dots represent 
the numerical solution, the dashed curve is the lowest order approximation given by Eq. (31), and 
the solid curve includes the next order approximation given in Eq. (32). 
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FIG. 11: As in Fig. (10), but for the reaction N + jcmb^^ ^ AT' + tt 

Similarly, [ebc~"'~^ f{a)] ^Z^""*"^-* is of order unity for a — ai* and rapidly increases for a > a*i. 
Thus for a sufficiently above oi*, one root (xi) is very close to unity and the other root (^2) 
is very large. Expanding Eq. (21) correspondingly, we find the following approximations for 

the two roots: 

1 - (a + l)e6c-°-V(a) 



Xi 



1 - (a + 2)e6c-«-V(a) ' 



(35) 



X2 ^ [ebc-^-^fia)] 



-l/(a+l) 



(36) 
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e < 



For e negative, Eq. (21) has one root, Xq. Using the same arguments as in the previous 
subsection, we conclude that for a sufficiently below a*i, the root xq is very small; and for a 
sufficiently above oi*, the root xq is very close to unity. Corresponding expansions give 



Xo 



[|e|6c-«-^ 

l+(a+l) |e| bc-^-^fja) 



/(a) J ' ; a < a*i , 



(37) 



a > a*i . 



l+(o+2) |e| 6c-°-i/(a) ' 

Figs. 12 and 13 compare the lowest order approximations [Eqs. (35), (36) and (37)] for 
the roots as a function of a with the corresponding numerical solutions, for the 7 and n 
reactions, respectively. 




FIG. 12: Threshold energy in logig (-E'th/eV) as a function of a, for the reaction 7 + 77^5— *-e+ + e~. 
The dots represent the numerical solution, and the curves are approximations given by Eqs. (35), 
(36), and (37) (xq above Ocrit and xi are indistinguishable on the plot). The curves are drawn 
slightly beyond their applicability ranges for easy identification. The parameter e is set to ±1. 



Stochastic Distributions 

The results of the previous subsections for fixed e can be folded into a distribution in e in 
a straightforward manner. When the variable e is distributed according to the probability 
density function p(e), the corresponding probability density function p{x) for the variable 
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FIG. 13: As in Fig. (12), but for the reaction N + 7cmb^A ^ N' + Tr 



de 



dx 



(38) 



X — x{€) is found from 

p(x) =p(e) 

Two of our four solutions, Xq with a > acrit in Eq- (37) and Xi in Eq. (35), are nearly 
independent of e. So for these two solutions, whether e is fixed or distributed is barely 
relevant. However, the other two solutions, xq with a < Ccrit in Eq. (37) and X2 in Eq. (36), 
are quite dependent on e. Taking a normalized Gaussian distribution 



p(e) = (27r(72)-2exp(-eV2(72) 



(39) 



as in the main text, and considering the roots xq for a < a*i and X2 for a > a*i, we find 



p{x) 



(27ra2)- 



(27r(7')- 



(« + 2)c"+i 
hf{a)x 



a+3 



exp 



. 1 (a + l)c«+i 



exp 



2(j262/2(a)x2«+4 

„2a+2 



2aWP{a)x 



2a+2 



, a < (3^,1 ; 



, a > a*i . 



(40) 



6/(a)a;«+2 

A few properties of the distribution p{x) are evident from Eq. (40): (i) p{x) is asymmetric 
since it decreases exponentially (as a power) to the left (right) of its peak; (ii) in the two 
intervals, a < a*i and a > a*i, the location and the height of the peak oip{x) increases with a; 
(iii) in both intervals, the width of the distribution decreases with a. These three trends are 
evident in the figures of this paper. To better compare p{x) with its numerical counterpart 
obtained in the paper, one would need to transform p[x) into p(log]^o ^) = (in 10) xp{x) and 
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adjust the normalization coefficient. We stop short of this. 
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